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Executive Summary 


Climate is, in many instances, the dominant variable controlling the 
storage of carbon in soils. It has proven difficult, however to determine how soil 
properties influenced by climate, such as soil temperature and soil moisture, 
actually operate to determine the rates of accumulation and decomposition of 
soil organic matter. Our approach has been to apply a relatively new tool, the 
comparison of 14 C in soil organic matter from pre- and post-bomb soils 
[Trumbore, in press (appended)], to quantify carbon turnover rates along 
climosequences. This report details the progress made toward this end by work 
under this contract. 

1) We have completed studies of carbon turnover rates from sites located 
along elevation gradients in Hawaii and the western Sierra Nevada mountains in 
California. Elevation gradients were chosen as the most effective way to study 
climate influences on soils as they are constant in the soil forming factors of 
parent material and topography, but differ in vegetation and climate. Two 
manuscripts from this phase of the work are attached (one in preparation, 
summarizing both transects, and one in press in Ecology) have resulted from this 
study. One of these manuscripts (attached: Climate controls on soil carbon 
turnover), summarizes both climosequences, and the implications for response of 
soil carbon storage or release accompanying a temperature change, and serves as 
a summary of this research and its implications for the role of soils in the global 
carbon cycle. 

2) We have completed a thorough study of the efficiency of different 
physically- and chemically-based fractionation methods at separating soil 
organic matter into pools which turn over on annual, decadal and centennial- 
millennial time scales. A preliminary manuscript summarizing this work (to be 
submitted to Radiocarbon) is also attached. 

3) A separate modeling study to understand the controls of the 14 C- 
content of soil carbon dioxide was undertaken by Yang Wang, a postdoctoral 
researcher on this project. Two papers published from her work are attached. 

In addition, three graduate students and one undergraduate student at UCI have 
received part-time support from this grant. 
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Wang, Y., R. Amundson and S. Trumbore, 1993. Processes controlling the 14C 
content of soil carbon dioxide: Model development. Chemical Geology 
107:225-226. 

Wang, Y., R. Amundson and S. Trumbore, A model for soil ^ CO 2 and its 

implications for using 14 C to date pedogenic carbonate (in press) Geochimica 
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ABSTRACT 

Climate acts as one of the major controls of the storage and turnover of 
carbon in soils [Jenny, 1949; Jenkinson et al, 1991; Schimel et al, in press]. We 
investigated the importance of climate in determining the dynamics of soil organic 
matter by measuring carbon and 14 C in organic matter from soils spanning the 
range of elevations in the western slope of the Sierra Nevada National Forest, 
near Fresno, California. These soils form a 'climosequence', a series of soils 
which are similar to one another in soil forming properties such as parent 
material, soil age, and topography (relief, slope, sun angle), but experience 
different climate [Jenny et al, 1949]. The amount of soil organic matter in detrital 
and soil A horizon layers increases with elevation along this transect, primarily 
reflecting increases in the amount of low density and hydrolysable organic matter 
fractions (fast cycling carbon pools). Turnover times of these fast-cycling pools, 
derived from modeling the increase in 14 C in SOM since the end of atmospheric 
weapons testing in the early 1960's, show a marked decrease with increased 
elevation (decreasing temperature). The observed relation between temperature 
and organic matter turnover for the Sierra transect soils is similar to that 
observed in Histosols from an elevation transect in Hawaii [Townsend et al., in 
press]. Extrapolating this relation to predict the transient response of fast-cycling 
soil organic matter to a global temperature decrease of 1°C, we estimate a 
potential increase in the storage of carbon in fast-cycling pools of due to slowed 
decomposition rate of up to 1 .3 PgC in one year. 
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INTRODUCTION 

Much of the uncertainty in recent assessments of the global carbon budget 
stems from the inability to measure directly the influence of the terrestrial carbon 
reservoirs. Soil organic matter (SOM) contains roughly two thirds of terrestrial 
carbon, two to three times more carbon than is present as CO 2 in the 
atmosphere [Schlesinger, 1977; Post et al., 1981; Eswaran et al, 1993]. One of 
the difficulties of understanding soil organic matter is that it consists of a complex 
mixture of compounds that accumulate and decompose at rates varying from 
months to millennia. The question of most importance to determining the role of 
soils in the global carbon cycle is how much of this mixture exchanges with 
atmospheric CO 2 on timescales of decades and less. Recent work using carbon 
isotopic approaches [O'Brien and Stout, 1978; Harkness et al., 1986; Balesdent 
et al, 1987; Scharpenseel et al., 1989; Trumbore et al, 1989; Vitorello et al., 

1989; Jenkinson et al; 1992; Trumbore, 1993; Harrison et al., 1993; Townsend et 
al., in press], observations of carbon inventory changes in soils following land 
cover changes [Parton et al., 1987; Schimel, 1986; Davidson and Ackerman, 
1993] and ecosystem modeling approaches [Jenkinson and Raynor, 1977, 

Parton et al., 1987; Van Bremen and Feijtel, 1990; McGuire et al., 1992; Potter et 
al., 1993; Schimel et al., 1994], has demonstrated that more than half of the soil 
organic matter in the upper 20-30cm of soils has turnover times of less than 1 00 
years. Thus it is important to investigate carefully how the abundance and 
turnover time of these rapidly cycling components of SOM are controlled by such 
factors as vegetation, climate, soil texture and landform characteristics. 

Although SOM is an inhomogeneous mixture with a continuum of turnover 
times ranging from days to millennia, it is more tractable to model the soil by 
lumping components in to pools which turnover on annual (active pool), decadal 
to centennial (slow pool) or centuries to millennial (passive pool) timescales 
[Parton et al., 1987; Schimel et al, 1994]. The problem of quantifying soil organic 
matter dynamics then becomes one of determining the sizes and turnover times 
of these three pools, and understanding the physical and chemical bases 
underlying these divisions. Fast pools (mostly active plus some slow) determine 
the magnitude of annual carbon losses by decomposition, and thus at steady 
state will equal annual carbon inputs by detritus [Jenkinson et al., 1991; Schimel 
et al., 1994], 14 C, the radioisotope of carbon (half-life = 5730 years), is useful in 
separating fast (active+slow) from passive carbon pools in soil organic matter. 
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The 14 C content of SOM collected before 1960 primarily reflects the influence of 
passive soil organic matter, in which carbon has resided in the soil long enough 
for significant radioactive decay. Passive pools may be isolated from bulk soil 
organic matter using physical and chemical fractionation procedures [Paul et at., 
1964; Campbell et al., 1967; Martel and Paul, 1974; Goh et al., 1976, 1977, 

1984; Scharpenseel et al., 1968a,b; 1984; Trumbore et al., 1989, 1990;Trumbore 
and Zheng, submitted]. The incorporation of 14 C produced in the early 1960's by 
atmospheric thermonuclear weapons testing ('bomb' 14 C) into soil organic matter 
during the past thirty years provides a direct measure of the amount of 
slow+active soil organic matter [Trumbore, 1993; Townsend et al, in press, 

O'Brien and Stout, 1984; Harrison et al., 1993]. A time-dependent modeling 
approach is used to link 14 C measurements of soil organic matter sampled prior 
to 1960 with that of contemporary samples. A large increase in 14 C content over 
the past thirty years indicates that significant portions of the soil organic matter 
are exchanging carbon with atmospheric CO 2 on decadal and shorter timescales. 
A comparison of 14 C in pre- and post-bomb soils also provides the best means to 
test the usefulness of various operationally defined separation procedures. 
METHODS 

Figure 1 shows the location of the Sierra Nevada climosequence, first 
described by Jenny et al. [1949]. Temperature and precipitation changes with 
elevation for the sites are summarized in Table 1, together with a description of 
vegetation and soil carbon inventories. Mean annual temperature decreases with 
elevation according to the atmospheric lapse rate. Precipitation increases rapidly 
to about 1000 m, then levels off. Most of the precipitation at higher elevations 
falls as snow. Vegetation also changes dramatically with elevation, from annual 
grasses in the Sierra foothills to pine forests at high elevation. Although we will 
summarize modeling results in terms of temperature controls on decomposition, it 
should be noted that other variables, such as the structure of vascular carbon 
inputs to the soil and precipitation also vary along the transect. 

Soils sampled in 1958-1962 as part of the California Vegetation Survey 
were available from the archive at UC Berkeley. The samples were stored air-dry 
in sealed glass Mason jars. No samples of detritus were archived. We used the 
original field notes and topographic map locations to find and re-sample the same 
sites in 1992. Additional data on litterfall and detrital mass measured along the 
same climosequence are reported in Jenny et al [1949], although no data were 
available for elevations below 4000 m. 
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Soil carbon and surface detrital inventories were determined from %C and 
bulk density measurements (corrected for gravel contents of up to 35% at high 
elevation in soil A horizons). Per cent carbon data were determined by CO 2 
evolution on in vacuo combustion at 900°C with cupric oxide wire. 14 C analyses 
were made by Accelerator Mass Spectrometry at the Center for AMS, Lawrence 
Livermore Laboratory, Livermore, California [Southon et al., 1992] from the 
purified CO 2 . Graphite targets for AMS measurement were prepared using the 
sealed-tube, zinc reduction method modified from Vogel [1992]. 14 C data are 
reported as A 14 C, the per mil difference in 14 C/ 12 C ratio between the sample and 
an absolute standard (oxalic acid decay-corrected to 1950; Stuiver and Polache, 
1977). 5 13 C values used to correct for mass-dependent fractionation effects were 
assumed to be -25%o for all samples. 

Soil organic matter was separated by density into fractions less and 
greater than 2.0 g/cc using sodium polytungstate [Trumbore, 1993; Trumbore 
and Zheng, in preparation]. The low density fraction consists of relatively 
undecomposed vascular plant matter and a small quantity of charcoal. The 
dense fraction (>2.0 g/cc) is predominantly mineral associated organic matter. 
Further treatment of the dense fractions (>2 g/cc) with a sequence of acid-base- 
acid (0.5N HCI; 0.1 N NaOH- Na; hot 6N HCI) isolated a residue which was 
always depleted in 14 C compared to the starting material. This residue was 
assumed to represent the A 14 C value and minimum carbon content of the 
passive pool. We report carbon inventory and 14 C data for low density (LF) 
fraction, the hydrolysable organic matter with density >2.0 g/cc (HYD; calculated 
by mass balancing C and 14 C between the untreated dense fraction and the 
residue), and the hydrolysis residue (RES) in Table 2. 

RESULTS 

Changes in carbon inventory and 14 C with elevation are summarized in 
Tablel (litter), Table 2 (soil A horizon) and Figures 2a, b. The total amount of 
carbon in both litter and mineral soil A horizons increases substantially with 
elevation. The exception to this is the Chiquito soil (highest elevation), which 
showed evidence of surface erosion (rills and collection of detritus upslope of 
fallen tree trunks). Figure 2a compares the carbon inventory in surface detritus 
measured in this study with the data reported by Jenny et al. [1949], 

The inventory of carbon in the soil A horizons (Figure 2b) equaled or was 
less than that in surface detritus at all elevations. Most of the increase in total 
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carbon inventory with elevation is due to increases in the low density(LF) and 
hydrolysable components (HYD). The amount of residual carbon after hydrolysis 
of the dense fraction in acids and bases (RES) is highest in the Musick soil, 
which also has the greatest clay content of any of the transect soils. 

One of the assumptions implicit in our comparison of 14 C in pre- and post- 
bomb soils is that all organic matter pools have remained at steady state over the 
past 30 years. The carbon inventory data given in Table 1 and Figure 2b are for 
1992 soils. Differences of carbon inventory between 1959 and 1992 profiles in 
mineral associated organic matter are less than 10% of the 1992 values for 3 of 
the 5 soils. The Chiquito, 1959 and Shaver, 1959 dense fractions had higher 
inventory than 1992 values by 100% and 65% respectively. Differences in LF 
inventory between 1959 and 1992 samples (which vary by less than 50% of the 
1992 values in all cases) may be due to differences in sampling and in the 
treatment of samples prior to storage (for example, by selective removal of larger 
roots). The fractional distribution between hydrolysable and nondydrolysable 
components was the same for both soils, though total inventory of dense material 
differed between 1959 and 1992 soils. We assume the carbon isotopes within 
each modeled component (LF, HYD, RES) will be independent of the relative 
amount of the component present, making spatial heterogeneity less important 
within a soil unit than the relative mix of LF, HYD and RES components. 

Detrital Layers 

Jenny et al [1949] report litterfall data based on one year of litter 
collections at the sites above 1200 m elevation along this transect. Turnover 
times calculated by dividing carbon inventory by annual litterfall [Jenny et al., 
1949], ranged from 8 to 115 years, and showed significantly faster turnover for 
deciduous than evergreen detritus (see Figure 3). As litterfall data were not 
available for the sites sampled in 1992, we assumed litter could be approximated 
as a homogeneous pool with a single average turnover time. We then compared 
the 14 C content observed in homogenized litter in 1992 to values predicted for a 
single pool at steady state with varying turnover times [see Harrison et al., 1993]. 
The measured 14 C data are compared in Figure 4a, with 1992 14 C values 
predicted using the turnover times derived for the detrital layers measured by 
Jenny et al. [1947]. Good agreement between the two approaches and data sets 
is observed at mid-elevations, while at high elevations, the detrital layer in the 
1992 Chiquito soil is probably affected by erosion (which causes an 
underestimate of turnover time). Turnover times for detritus based on 14 C will 
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include the residence time of needles and branches on trees (which may be up to 
several years for evergreens at high elevations). 

Soil A Horizons 

Figures 4b, 4c, and 4d show the 14 C increase between 1959 and 1992 for 
A horizon LF, HYD and RES fractions. Fractions with slowest turnover are 
expected to have the least 14 C initially, and the smallest increase in 14 C over the 
past 30 years. Based on these criteria, the RES fraction has the longest turnover 
time (and will be used in modeling as representing the passive pool). Low density 
carbon (LF; Figure 2b) shows the largest increase between 1959 and 1992 (the 
Fallbrook soil was collected in 1962. The 14 C values for atmospheric CO 2 in the 
1950's and in 1992 (and therefore plant carbon fixed during these years) are 
indicated on Fig 4b by arrows at roughly -25%o and +140%o [Stuiver and Quay, 
19xx; Camargo, in preparation], 14 C values of LF in archived soils were less 
than the 14 C content of the atmosphere in the 1950's (which shows the effect of 
dilution with 14 C-free fossil-fuel derived CO 2 ; Suess, ), indicating the presence of 
some inert or passive fraction (most probably charcoal, see below). 14 C values in 
1992 LF generally decrease with elevation, ranging from slightly above present- 
day atmospheric values in the Fallbrook soil (lowest elevation) to below-1992 
values at higher elevations. Fast turnover of LF at the low elevation site 
(Fallbrook) is evidenced by the presence of significant amounts of bomb 14 C in 
the 1962 archived soil. 14 C values in the hydrolysable dense organic matter 
(HYD fig 4c) are intermediate between RES and LF fractions. 

Models 

A three-component model was used to determine the distribution and 
turnover rates of carbon in active, slow and passive components needed to 
reproduce the observed increase in 14 C for density <2.0 and >2.0 fractions 
[described in detail in Townsend et al., in press; Trumbore et al., in preparation]. 
Turnover times in the model were constrained to be one year for the active pool, 
and by the observed 14 C content of the RES fraction in archived soils for the 
passive pool. The amount of RES carbon was also used to set the minimum size 
of the passive pool. The turnover time of the slow pool was then adjusted until 
both archived and 1992 14 C values were matched. An additional constraint 
(which limited the amount of active SOM to close to zero everywhere but in the 
Fallbrook soil) was the total annual carbon flux through the soil. One estimate for 
this is obtained from the above-ground litterfall data of Jenny et al [1949], which 
ranged from 80 to 160 gC nrr 2 a* 1 (1300 to 3000m elevation). A second estimate 
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was derived using the relation of temperature, precipitation and CO 2 flux from 
Raich and Schlesinger (1991). We assumed 50% of the total soil respiration, or 
roughly 200-300 gC nr 2 a -1 as calculated using Raich and Schlesinger (1991) for 
these soils, was due to CO 2 released by decomposition. 

We modeled LF (density <2.0) and HF (density >2.0 g/cc) fractions 
separately, as each contained slow and passive components. We derived total 
turnover times for active, slow and passive pools by adding together the amount 
of active, slow or passive carbon in LF with that in HF, and dividing the sum for 
each pool by the sum of annual carbon flux through it (calculated by summing 
contributions to the total flux (amount/turnover time) from active, slow and 
passive components of LF and HF). As only the Fallbrook soil showed a 
significant amount of carbon in the active pool, and to facilitate comparison with 
reults from other locations and models (see below), active+slow pool turnover 
times were combined in the same matter. The results are summarized in Table 
3. 

The modeling of the LF fraction was complicated by the presence of 
charcoal. As there were no fires at these sites between 1959 and 1992 
collections, any charcoal present in the soil will contain pre-bomb carbon, which 
serves to dilute 14 C content of vascular plant remain components of LF. Visual 
inspection of the LF material revealed charcoal, but no attempt was made to 
quantify what proportion of the total carbon was in this form. In order to match 
both 1959 and 1992 LF 14 C values, between 5 and 15% of the carbon was 
required to be inert, with 14 C values of -50%o. Large amounts of active carbon 
(constrained as having one year turnover) were precluded by the high inventories 
of carbon in low density fraction, coupled with the knowledge that carbon inputs 
are probably of the order of 5% or less of LF carbon inventories. 

The dense fraction 14 C content could be modeled reasonably well using 
the abundance and 14 C content of the RES fraction as representing the passive 
pool. The HYD fraction was then nearly all in the slow pool. Turnover rates in the 
HF slow pool were equal to or slightly greater than those in the LF slow pool for 
the same soil. 

The turnover times for combined active+slow A horizons are plotted with 
surface detrital values in Figure 3. Fast cycling components of A horizon soil 
organic matter show an increase in turnover time with elevation similar to that 
observed in surface detritus. 
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The utility of the physical and chemical fractionation procedures at 
separation of SOM directly into active, slow and passive constituents is assessed 
in Table 3. The sum of active +slow carbon inventory is 77 to 98% of the sum of 
LF + HYD. Thus LF+HYD is a good predictor of the amount of active+slow 
carbon in these soils. This is not true for other soils, such as Histosols, which do 
not show large differences in 14 C on fractionation (Goh et al., 1977, 1984; 
Trumbore and Zheng, submitted). 

Discussion 

Figure 5a plots the turnover times of the active plus slow pools from the 
Sierra Nevada climosequence against mean annual temperature. The results 
from another elevation-based climosequence on Histosols in Hawaii [Townsend 
et al, in press] are also shown, along with a datum representing an Oxisol from 
the Amazon Basin, Brazil [Trumbore et al., in preparation]. Although these soils 
differ widely in vegetation cover, precipitation, and soil parent material, the strong 
relation between turnover time of the combined active + slow pools and 
temperature demonstrates that it is strong controlling factor on the dynamics of 
carbon in soils, as argued in recent modeling and data-based studies [Jenkinson 
et al; 1991, McGuire et al., 1 992; Buol et al., 1 993; Potter et al., 1 993; Schimel et 
al., 1994], Litter turnover times derived from 14 C for the Sierra soils also show a 
relation of decomposition to temperature similar to (or steeper than) that for soil A 
horizon active+slow pools (see Figure 3). 

The carbon inventory contained in active+slow pools in these soils is 
plotted in Figure 5b. At each site, the active +slow carbon makes up more than 
60% of the total carbon in the soil A horizon, and carbon inventory within each 
climosequence is inversely related to temperature. The slopes derived from the 
two climosequences are similar (excluding the Chiquito soil in the Sierra 
transect), although the absolute amount of carbon is much greater in the Hawaii 
climosequence soils. 

The temperature relation to carbon abundance and turnover in active + 
slow pools observed in these regional scale elevation gradient studies are 
compared in Figures 5a and 5b with results derived from global ecosystem 
models [Potter et al., 1994; Schimel et al; 1994], The carbon inventory and 
turnover rates of active+slow soil carbon pools derived from the CASA model 
[Potter et al., 1994, and Potter, personal communication, 1994], averaged into 
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10° latitude bands, are plotted using small cross symbols. The model-derived 
turnover times are a less strong function of mean annual temperature than the 
climosequence-based relation. Exponential fits to the climosequence (+Brazil 
soil) data give the following relation between turnover time and man annual 
temperature (MAT): 

Turnover(A+S) = 74*exp(-0.074*MAT) 1), 

while a fit to the model-generated points is: 

Turnover(A+S) = 41*exp(-0.0546*MAT) 2). 

The importance of soil organic matter and detrital pools in sequestering or 
releasing carbon given a temperature change may be estimated using equations 
(1) and (2), and the estimates of carbon storage by latitude taken from Potter et 
al [1994], We assume that a transient temperature change, such as a decrease 
in MAT caused by a large volcanic eruption, will act to slow carbon losses from 
soil by decomposition, while leaving carbon inputs unchanged. The transient 
(one year of forcing) response of the global active+slow carbon pools is then 
estimated using the data in Table 4 for carbon inventory, and using turnover 
times calculated using equations 1 and 2: 

A(carbon storage) = Inputs - Outputs 

= C inventory/turnover(MAT) - C inventory/turnover(MAT-1) 

The data used for this calculation (Table 4) are lumped into 5 latitude bands, with 
distribution of land area by latitude derived from Potter et al. [1994]. 

The results for both the climosequence and model-derived relations 
for turnover vs. MAT show the potential for one-year storage of 1 .2 to 1 .3 Pg 
of carbon in active+slow pools of soil A horizons with a 1°C decrease in 
MAT. The largest changes in soil carbon storage occur at low latitudes (high 
MAT). Change in the carbon inventory in any given soil would be 
undetectable (<3% of A horizon inventory). Further storage would 
presumably occur in surface detrital layers, which have a similar or steeper 
response of turnover rates to MAT (Fig 3). These calculations ignore some 
important complicating factors, such as a probable slowing of carbon inputs 
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to the soil by vegetation at lowered MAT, and feedbacks between slowed 
decomposition and plant productivity (Schimel et al., 1994). In addition, the 
climosequences studied integrate long-term adjustment of ecosystems to 
climatic conditions, and transient responses may not be predictable from 
climate sensitivity derived from these data. However, they serve to illustrate 
that changes in MAT can potentially cause a globally important transient 
response by storing or releasing carbon from active+slow carbon pools. This 
response may provide a partial explanation for a slowing of the rate of 
atmospheric CO 2 increase in 1993, associated with global cooling following 
the eruption of Mt. Pinatubo [NOAA Climate Atlas, 1993]. 

Conclusions 

Temperature is major control on the dynamics of organic matter which 
turns over on decadal and shorter timescales in soil A horizons. The relationship 
derived for turnover time vs temperature from climosequence studies is more 
sensitive to mean annual temperature changes than a relationship derived from 
the CASA global ecosystem model [Potter et al., 1994, Schimel et al; 1994], 

The data presented here represent only two climosequences, from soils 
with very different vegetation, precipitation and temperature. Although more data 
are needed to provide a better understanding of climate controls on carbon 
turnover in soils, the data presented here demonstrate the usefulness of the 
climosequence approach. In particular, a better understanding of controls on 
carbon turnover in tropical soils is needed, as the response there is most 
sensitive. Comparison of the highest temperature points from climosequences in 
Figure 5a suggests considerable textural or vegetation controls may be important 
in these soils. 

Fractionation of soil organic matter, while not a perfect indicator of the 
distribution of turnover times, is useful indicator of which physical portions of the 
soil organic matter are more labile and more refractory. In the Sierra transect 
sequence, low density and hydrolysable mineral-associated organic matter made 
up the active and slow pools (with the exception of charcoal present in the low 
density fraction), while the residue left after hydrolysis of mineral-associated 
organic matter in acids and bases provided a reasonable estimate for the size 
and turnover time of passive soil organic matter (once again, with the exception 
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of charcoal). These fractionation procedures appear to be most useful in soils 
with sandy or silty parent material; clay dominated soils, such as the Hawaii 
climosequence soils, do not show large differences in 14 C among these particular 
fractionation methods [Townsend et al., in press; Trumbore and Zheng, in 
preparation]. 
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Figure Captions 

Figure 1. Location of Sierra Nevada Transect in Fresno County, CA. Approximate 
locations of the 5 soils on the elevation transect are also shown. 

Figure 2. Carbon inventory in surface detrital (litter) layers (Fig. 2a), from this 
study and from Jenny et a [1947], Figure 2b. Carbon inventory changes with 
elevation for the bulk A horizon soil organic matter (total), fractions with density 
<2.0 (LF), and hydrolysable (HYD) and nonhydrolysable components of the 
fraction with density >2.0 g/cc. 

Figure 3. Changes in detrital layer turnover time with elevation along the Sierra 
transect. Turnover times from data in Jenny et al [1947] (squares) were 
calculated by dividing the detrital inventory by the annual litterfall (both were 
assumed to have the same %carbon). Turnover times for 1992 Sierra transect 
litter were determined from the 14C content of homogenized litter, assuming that 
litter may be represented using a single turnover time and is at steady state (Fig 
4a). Estimated turnover times for A horizon active+slow components (derived 
later) are shown for comparison. 

Figure 4. Comparison of carbon-14 contents for organic matter fractions isolated 
from Sierra Nevada transect soils (1950's and 1992). Negative A 14 C values 
indicate that carbon in a given fraction has had sufficient time to undergo 
radioactive decay (half-life of radiocarbon is 5730 years), while positive values 
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indicate the presence of bomb-produced 14 C. Lowest initial 14 C values and the 
smallest increase in 14C over the past 30 years are observed in the RES 
fraction; highest initial 14C and greatest 14 C increase occur in the LF fraction. 

Figure 5a. Estimated turnover times in active+slow pools (see Table 3) from 
modeling of 14 C increases in density <2.0 (LF) and density >2.0 fractions. Data 
from a similar elevation transect study in Hawaii [Townsend, Ph.D. dissertation, 
and Townsend et al., 1994 in press], and a single datum representing carbon 
dynamics derived from 14C and carbon flux measurements in the eastern 
Amazon basin [Trumbore et al., in preparation; Camargo et al., in preparation]. 
Crosses (identified in the legend as GLOBAL) are derived from results of the 
CASA global ecological model [Potter et al., 1994; and Potter, pers. comm.]. 
Curve fits used for equations 2 and 3 in the text are also shown. Figure 5b. 
Carbon inventory in active+slow pools for the same soils and model results. 
These represent >50% of the total soil carbon in all cases. 



Table 1. Climate and vegetation of Sierra Nevada study sites 
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Table 2. Carbon inventories and 14 C in fractonated soil organic matter. 
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Table 3. Model results partitioning carbon into Active, Slow and Passive pools. 




Table 4. Predicted Chnage in Carbon Inventory from a 1°C change 
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Comparison of Fractionation Methods for Soil Organic Matter 14 C Analysis 
Susan E. Trumbore and Zheng Shuhui 

Department of Earth System Science, University of California, Irivne, CA 92717 

(prepared for submission to Radiocarbon) 

ABSTRACT 14 C measurements provide a useful test for determining the degree 
to which chemical and physical fractionation of soil organic matter are successful 
in separating labile and refractory organic matter components. Results from 
AMS measurements of fractionated soil organic matter made as part of several 
projects are summarized here, together with suggestions for standardization of 
fractionation procedures. While no single fractionation method will 
unequivocally separate soil organic matter into components cycling on annual, 
decadal and millennial timescales, a combination of physical (density separation 
or sieving) and chemical separation methods (combined acid and base 
hydrolysis) provides useful constraints for models of soil carbon dynamics in 
several soil types. 

INTRODUCTION 

Models describing the dynamics of accumulation and turnover of organic 
carbon generally recognize components of soil organic matter (SOM) which turn 
over on annual (active), decadal (slow), and centennial to millennial (passive ) 
time scales [Jenkinson and Raynor, 1977, Parton et al., 1987; Van Bremen and 
Feijtel, 1990; O'Brien and Stout, 1978; McGuire et al., 1992; Potter et al., 1993; 
Schimel et al., 1994]. While these concepts have proven useful in explaining the 
magnitude and timing of changes in SOM following a perturbation such as land 
clearing for cultivation [Parton et al., 1987; Schimel, 1986; Cambardella and 
Elliott, ;Davidson and Ackerman, 1993], no recognized method now exists for 
determining how to apportion SOM into compartments which turn over on 
different time scales. To be useful for dividing soil organic matter into 
components useful for ecosystem modeling, a procedure must determine both 
the pool sizes and turnover rates. 

Two methods are presently in use to partition soil organic matter into 
active, slow and passive pools. The first models the observed increase of 14 C 
during the 30 years since the end of atmospheric weapons testing, using 
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additional constraints derived from knowledge of carbon inputs to the system 
[O'Brien and Stout, 1978; Harkness et al., 1986; Balesdent et al, 1987; Scharpenseel 
et al., 1989; Trumbore et al, 1989; Vitorello et al., 1989; Jenkinson et al; 1992; 
Trumbore, 1993; Harrison et al., 1993; Townsend et al., in press]. The second 
relies on physical and/or chemical fractionation methods to separate SOM into 
pools which turn over on different time scales [Paul et al., 1964; Campbell et al., 
1967; Martel and Paul, 1974; Goh et al., 1976, 1977, 1984; Anderson and Paul, 

1984; Scharpenseel et al., 1968a,b; 1984a,b; Trumbore et al., 1989, 1990; Trumbore 
et al., in preparation]. Radiocarbon measurements of pre-bomb SOM fractions 
are interpreted as indicative of their turnover rates. A test of the efficiency of any 
fractionation method at isolating a relatively homogeneous (in terms of turnover) 
organic matter pool is a comparison of the observed 14 C increase since 
atmospheric weapons testing (Goh et al., 1984; Trumbore et al., 1989; Trumbore 
et al., in preparation). These tests have shown that, while no fractionation 
method is completely successful, separation of OM by density, followed by 
hydrolysis in 6N HC1, does result in a collection of compounds which are on 
average, turn over more slowly than the hydrolyzed portions of SOM [Campbell 
et al., 1967; Martel and Paul, 1974; Goh et al., 1976, 1977, 1984; Anderson and 
Paul, 1984; Scharpenseel et al., 1968a, b; 1984a,b; Trumbore et al., 1989, 1990]. 

This paper will summarize measurements of * 4 C made of soil organic 
matter fractionated using different physical and chemical separation methods. 
Both archived (pre-bomb) and contemporary soil samples have been analyzed. 
The implications of these data for soil carbon cycling is discussed in other 
publications (see Table 1 for references). This paper is in essence a 'date list' 
which will make the original data available to those interested in the effects of 
fractionation methods on the 14 C content of SOM. 

Summary of Fractionation Methods 

All samples are sieved to remove components >2mm. Chemical and 
density separation procedures have been described before (Trumbore et al., 1989; 
Trumbore, 1993), and are summarized in Figure 1. Basically three kinds of 
fractionation techniques are used: separation by density, chemical extraction 
using acids and bases, and separation by size. 

1. Density separation. The density separation extracts organic matter of low 
density (<1. 6-2.0 g cm-3) by flotation in a heavy liquid. The low density material 
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consists of plant matter and charcoal [Sollins et al, 1983; Spycher et al, 1983]. 
Dense fractions contain mineral-associated organic matter, and some microbial 
cell debris [Sollins et al, 1983; Spycher et al, 1983]. In general, the low density 
material turns over more rapidly than much of the soil organic matter. 

2. Chemical Fractionation. Operationally defined fractionation procedures for 
SOM rely on extraction using acids and bases, which are hypothesized to remove 
more labile components, leaving behind more refractory (and therefore more 
14 C-depleted) constituents. The data presented here compare the 14 C content of 
residues after hydrolysis in acid (HC1) of varying strength, and in base (NaOH or 
NaOH-Na-pyrophosphate). Because the extraction of fresh plant material by 
these techniques leaves a residue [Ertel and Hedges, 1984], we have chosen to 
perform acid and base hydrolysis on SOM which has been separated previously 
by density, and to treat only the denser portion of the soil further. 

3. Size separation. Recent work has tied the dynamics of decadal-cycling SOM 
to different size fractions (Cambardella and Elliott, 1993; 1994), which has led to 
speculation about the role of soil aggregates in limiting the decomposition of 
reactive SOM constituents by physically protecting organic matter from soil 
microbial activity [Oades, 1993]. In addition, the limited measurements of 14 C in 
size fractionated organic matter available from the literature [Anderson and Paul, 
1984] suggest a relationship between the size of soil-associated organic material 
and stability of associated organic matter. Recently, Mayer [1994] has shown that 
organic matter content is related to BET surface area in many soils and 
sediments. If organic matter is protected by stabilization on clay surfaces, we 
might expect it to have lower 14 C values than organic matter which is more 
available for microbial attack. 

For the size separation, we used the <2mm sieved samples, without prior 
density separation. We measured the 14 C content of three size fractions: <2mm 
but >63|i, <63 (i but greater than 2p, and less than 2p. The largest fraction was 
separated by dry sieving; the fractions > and < 2 p were separated by settling 
from a suspension [Folk]. This procedure breaks apart soil aggregates. 

Per cent carbon data were determined by CO 2 evolution on in vacuo 
combustion at 900°C with cupric oxide wire [Buchanon and Corcoran, 1959]. 14 C 
analyses were made by Accelerator Mass Spectrometry at the Center for AMS, 
Lawrence Livermore Laboratory, Livermore, California [Southon et al., 1992] 
from the purified C02- Graphite targets for AMS measurement were prepared 
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using the sealed-tube, zinc reduction method modified from Vogel [1992]. 14 C 
data are reported as A 14 C, the per mil difference in 14 C/ 12 C ratio between the 
sample and an absolute standard (oxalic acid decay-corrected to 1950; Stuiver 
and Polach, 1977). 5 l3 C values used to correct for mass-dependent fractionation 
effects were assumed to be -25 %o for all samples. 

Nomenclature: We adopt a shorthand for describing the fractionation treatment 
experienced by a sample. The descriptors, denoted in Figure 1, are defined as 
follows: 

BS = Bulk soil. Due the expense of AMS measurements, we often rely on 
measurements of the constituent fractions to derive a bulk soil 14 C number by 
mass balance. 

LD = low density fraction (less than 2.0 g/cc). 

HD = Dense fraction. The density of heavy liquid used in the separation is 
identified by accompanying LD data. 

Chemical fractionation. All chemical fractionation procedures are performed on 
the HD fraction, and the data associated with them refers to the residue after a 
given chemical separation procedure. A refers to acid hydrolysis at room 
temperature in 0.5 N HC1; B refers to extraction with 0.1 N NaOH- 0.1 N sodium 
pyrophospate; C refers to hydrolysis in 6N HC1 at 95°C. Often, we show data 
derived from the residue of a combination of chemical fractionation procedures: 
for example A+B+C refers to the residue after sequential hydrolysis in 0.5N acid, 
0.1N NaOH and 6N HC1. When performing these procedures, the residue after 
each extraction step is rinsed thoroughly in distilled H 2 O before moving on to 
the next procedure. 

Samples 

We present here summaries of fractionation of soils from several different 
types. Table 1 summarizes the source of soils, and data related to soil 
classification and site-specific characteristics, such as climate and parent material. 
References discussing the carbon cycling in these soils, or containing more 
information about the soils, are also listed. 


Results 
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Tables 2 and 3 show the carbon and 14 C content of the fractionated soil 
organic matter. To better compare the effects of fractionation methods on the 
different soils, the results are also plotted in Figures 2 and 3. 

Density separation 

Low density components consist primarily of fine root hairs, charcoal 
(found in all LD fractions, but not quantified separately), and surface detritus 
mixed into the soil (in Oxisols BS-7, BS-9). Low density material makes up a 
large portion of the total soil carbon in A horizons of more coarse-textured soils 
and the Histosol (which had the highest total carbon content. Table 2). The 14 C 
content of LD material is close to atmospheric 14 CC >2 values at the time of 
collection in both pre- and post-bomb soils, indicating rapid turnover of this 
fraction. 

In general, 14 C values of the dense fraction (mineral associated) carbon 
decreased with depth in the soil, as has been previously observed (Scharpenseel 
et al., 1989). Bomb 14C is present in HD fraction for soils sampled in the 1990's, 
indicating the presence of at least some fast-cycling carbon associated with soil 
mineral phases. 

Acid-Base hydrolysis 

In the more heavy-textured soils (those excluding the Shaver soil from the 
Sierra Nevada, California), most of the carbon present is associated with mineral 
surfaces (HD fraction). Hydrolysis of the HD fraction in relatively mild acids 
(0.5N HC1) removes greater than 30% of the mineral-associated carbon for all 
soils studied (see Figure 2a). The removal is greatest in Spodic B horizons (the 
Michigan soils samples, Bl-56 and Al-52), where more than 70% of the organic 
matter is solubilized in 0.5 N HC1. The carbon remaining is depleted in 14 C 
compared to the starting dense fraction material to varying degrees (Fig 3a). 

Least affected are the Histosol (Hawaii) samples, most affected are the Spodic B 
horizons. Subsequent hydrolysis in base and strong acid removes more carbon 
and 14 C (Fig 2b and Fig. 3 a-c). The residue after A+B+C treatments is, in all but 
two cases, the most depleted in 14 C of any treatment. Figure 3 summarizes the 
difference in 14 C content of residual carbon after various treatments, compared to 
the starting dense material (1:1 line on Figures). Use of A+B treatment (Fig 3b) 
removes almost as much carbon and 14 C as the A+B+C treatment, while 
treatments using only acids (A+C; Fig 3c) do not remove as much C and 14C as 
treatments including a base hydrolysis step. The Histosols and Oxisols are 
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affected by chemical fractionation less than more coarse textured soils. The 
difference between unfractionated dense material and the A+B or A+B+C 
residues is smaller in A horizons of soils than in B horizons. This is more 
noticeable in profiles previously published in Trumbore et al. [1990]. 

Size fractionation 

The results of size fractionation on splits of the same soils in which 
density- and chemical procedures were used are reported in Table 3. Most of the 
carbon in more coarse textured soils was associated with the > 63|i fraction 
(Figure 2b), while in heavy-textured soils (Oxisol), most of the carbon is in the 
fine silt and clay fractions (<63 ji). Size separations were not performed on the 
Histosol (Hawaii ) soils; these are dominated by fine-grained material. Size- 
fractionation of soils did not affect 14 C values as much as chemical fractionation 
(Fig 3d). We did not see any consistent tendency for clay-sized particles to have 
A 14 C values much greater or less than other size fractions, although clay fractions 
both enriched and depleted in 14 C were observed. 

Conclusions and Suggestions for a standard fractionation scheme. 

The usefulness of different fractionation methods depends in part on the 
type of soil, in particular on soil texture (clay content). Separation by density will 
concentrate organic matter which turns over fast (decades and shorter timescales) 
in the low density fraction. In coarse textured soils with granitic parent material, 
such as the Shaver soil (Sierra Nevada), low density organic matter can make up 
the majority of the total soil carbon. The remaining HD material shows 
significant increase in 14 C in A horizons of contemporary soils due to 
incorporation of bomb 14 C. Chemical treatment, especially hydrolysis in acid 
and base, removes relatively 14 C-enriched components of the dense, mineral- 
associated organic matter. Once again, chemical treatment is most effective in 
coarser textured soils, in part because more of the total carbon is removed in the 
hydrolysis procedure. Chemical extraction with acids and bases did not change 
the 14 C content of Histosol organic matter, even though up to 90- 95% of the 
original carbon was removed. This may be useful in understanding how organic 
carbon is stabilized in allophane soils. 

Large differences in ^ 4 C content were not observed between different size 
fractions. Most of the carbon in coarse-textured soils is in the larger size fractions 
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(>63p), and probably is mostly low density vascular plant material, though 14 C 
values indicate that this organic matter is diluted by other (mineral associated) 
constituents. 

We suggest the most informative fractionation procedure to be used in 
separating soil organic matter into faster and slower cycling components 
combines density and chemical separation procedures. Low density organic 
material (except for charcoal) turns over on timescales of decades and less (see 
Trumbore, 1993; Trumbore et al., in preparation; Camargo et al., in preparation). 
Extraction of dense material with acid and base will leave a residue which can 
provide a minimum estimate of the size and 14C content of a passive pool in 14 C 
modeling. This fractionation will be most useful in soils of medium to coarse 
texture. 
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FIGURE CAPTIONS 

Figure 1. Fractionation scheme used for 14C measurements. 

Figure 2. Per cent of total carbon in the residues left after acid and base 
hydrolysis treatments of the dense fraction (Fig 2a), and after size fractionation of 
bulk soil (Fig 2b). 

Figure 3. A 14 C of HF (x axis) compared to that measured in the residue after 
various chemical treatments (identitifed in panels). The 1:1 line and the A14C 
values for A+B+C (most 14 C-depleted treatment) are given on each figure for 
comparison. 



TABLE 1. Summary of Characteristics of Soils Analyzed in the Fractionation Study. 
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Table 2 . Results of Chemical Fractionation Procedures 
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assumed 5 13 C of -25%o for the purpose of calculating A 14 C for all samples, all vegetation cover for these soils was C 3 plants. 
The one sigma error of the 14 C measurements is ± If) % c to 15 %c. 



Table 3. Results of Size Fractionation Procedures 
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A model for soil l4 C0 2 and its implications for using ,4 C to date pedogenic carbonate 
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Abstract-A diffusion-reaction model for so.l l4 C0 2 is described that analyzes the “(^distribution in 
S0 ‘ ^' t sh h °^ s tha . t . the C content of soil C0 2 is not the same as that of atmospheric C0 2 and varies 
with depth depending on various factors. The most important factors affecting the 14 C content of soil 
rvV mc ' ude lhe content of soil organic matter, the relative contribution of root respired CCU to total 
C0 2 production, soil respiration rate, atmospheric C0 2 concentration and l4 CO, content and soil properties 
such as temperature and moisture content etc. The l4 C content of soil CO, not onlv can be a sensitive 
indicator of the residence time of decomposing organic matter in the soil, but also determines the I4 C 
content of pedogenic carbonate. Our model suggests that soil C0 2 could be enriched or depleted in l4 C 
relative to atmospheric C0 2 , depending on the relative contribution of root respiration to total soil 
respiration and on the turnover rate of the soil organic matter contributing to the soil C0 2 . Therefore, 
the initial C dates of soil carbonate could differ from the true ages of pedogenesis. The processes and 
factors considered by the model are a first step in determining whether the 14 C content of soil carbonate 
could lead to reliable dates of pedogenesis. 


INTRODUCTION 

Carbon dioxide IS PRODUCED by decomposition of organic 
matter and by root respiration in soils. As a result, CO : con- 
centrations in soils are higher than atmospheric CO : con- 
centrations. Diffusion is considered the major mechanism by 
which gases move in soils (Kjrkham and Powers, 1972: 
Jury et ah, 1991). It has been previously shown that the 
distribution of CO : and its stable isotopic species (e.g., l2 CO : > 
13 C0 2 ) in soils can be described by a diffusion-reaction model 
(Cerling, 1984, 199 1 : Solomon and Cerling, 1987; Cer- 
LlNCetai., 1989; QuadE et al., 1989). Since the<5 l3 C value 
of C0 2 produced by root respiration is essentially the same 
as that of C0 2 produced by decomposition of soil organic 
matter, the S ,3 C value of soil C0 2 is determined by soil res- 
piration rate and the 6 l3 C value of soil organic matter, w hich 
is more or less constant with depth and time in a stable soil 
ecosystem (Cerling, 1984, 1991; Cerling et al., 1989; 
Quade etal., 1989). 

However, in the case of the radioactive carbon isotope 
( M C), the model needed to describe the l4 C0 2 profile will be 
more complicated. First of ail, the 14 C content of CO^ pro- 
duced by root respiration is not necessarily the same as that 
of C0 2 derived from organic matter decomposition. While 
root respiration C0 2 can be assumed to be nearly identical 
with the atmospheric I4 C level (DOrr and MOnnich, 1986), 
the ,4 C content of C0 2 derived from organic matter decom- 
position may, depending on the carbon residence time in the 
soil, be depleted by radioactive decay or even enriched by 
the contribution of ,4 C produced in the late 1950s and early 
1960s by atmospheric testing of thermonuclear weapons 
(bomb ,4 C) (DOrr and MOnnich, 1986). In other words, 
the S !4 C value of CO : produced by organic matter decom- 
position should be determined by the 14 C content of soil or- 
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ganic matter, which is a mixture of compounds with different 
turnover rates (Trumbore et al.. 1990: Trumbore, 1993). 
Moreover, the o ,4 C value of organic matter changes with 
depth and time (Goh et al.. 1976: O’Brien and Stout. 1978: 
Trumbore. 1993). suggesting that the CO : produced at dif- 
ferent depths should have different l4 C contents depending 
on the l4 C content of soil organic matter at that depth and 
also on the relative contribution of root respiration to the 
total C0 2 production. Haas et al. ( 1 983 ) and Thorstenson 
et al. ( 1983) observed lower 14 C contents in respired CO 
during winter time due to oxidation of older soil organic 
matter in the absence of significant root respiration. However, 
their model failed to describe the 14 CO : -depth profile observed 
in their study. DOrr and MOnnich (1986) also observed 
an annual l4 C variation of soil-respired CO : , which they sug- 
gested is controlled by the seasonally vary ing contribution of 
root respiration and C0 2 produced by microbial decompo- 
sition of organic matter. 

Presently, there is no adequate model describing the I4 C0 2 
distribution in soils. Previous studies on soil CO : , and 5 13 C 
values of soil C0 2 , and pedogenic carbonates indicate that 
the C0 2 and carbonate in a soil system are in isotopic equi- 
librium. This implies that any detrital carbonate dissolved in 
the soil is overwhelmed by soil C0 2 and the isotopic com- 
position of pedogenic carbonate is determined by the isotopic 
composition of soil C0 2 (Cerling, 1984, 1991; Cerling et 
al., 1989; Quade et ah, 1989; Amundson et al., 1989; Cer- 
ling and Quade, 1992). If this is the case, the 14 C content 
of pedogenic carbonate should also be determined by l4 C 
content of soil C0 2 and inherited dead carbon should not 
affect the age of soil pedogenic carbonate. Therefore, a better 
understanding of ,4 C0 2 distribution in soils would have very 
important implications not only in the study of organic car- 
bon cycling in soils, but also in the application of ,4 C to date 
soil carbonate. Furthermore, the l4 C content of soil CO, is 
an important parameter in modeling variations of atmo- 
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where 


Tabic 1 Parameters for the model soil desenbed^in the text 
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According to Cerling (1984), the steady-state solutions 
for diffusion-reaction equations for C0 2 and l3 C0 2 , with 
boundary conditions of C* = C* and Cl } = Cl 3 at z = 0, 
J dC* n dCl 3 

and — — — 0 and — — = 0 at r = L, are as follows: 
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The diffusion coefficients of 14 C0 2 and ,2 C0 2 are related 
by 




Dl 1 

7.W(air) + .V/( l2 C0 2 )\ 

/ .\f( air) X .V/( u C0 2 )\'| 

t) +cj 

(7) 

D\* 

\ .V/( air ) X .l/( l2 C0 2 )j 

\.V/(air) + .V/(' 4 C0 2 ) j 



where CJ and C[ 3 are CO : and ,3 C0 2 concentrations in the 
atmosphere: C? and C\ 3 are corresponding concentrations 
in the soil air; and D* and Z), 3 are diffusion coefficients for 
C0 2 and n C0 2 , respectively. 

The diffusion coefficient for C0 2 in soil is related to that 
in air <Z> a(r ) by 

D* = 

where t is the free air porosity in the soil, and p is a tortuosity 
factor (Kjrkham and Powers, 1972; Jury et aL 1991 ). 
£> air varies with temperature and pressure ( Bird et al.. 1 960 ) : 



where D * r is the diffusion coefficient for C0 2 in air under 
standard conditions ( T° = 25°C and P° = 1 bar pressure) 
and is taken to be 0. 144 enr/s. 


8 14 C (%c) 



= 1.00868. 

Using Eqns. 6, 7. and 8 and using various values of soil 
respiration, an atmospheric CO : concentration of 300 ppm. 
and atmospheric 5 l3 C and 6 ,4 C values of -6%o ( pre- industrial 
value ) and 0 %o ( pre-bomb value ). respectively, it is possible 
to calculate I4 CO : (6 M C) profiles in soils for various condi- 
tions. 

Figure l shows how the <5 U C value of soil CO, vanes with 
depth in a model soil where the I4 C0 2 concentration is dif- 
fusion-controlled. The parameters for the model soil are listed 
in Table 1. This diagram displays several important features: 

1 ) The 6 ,4 C value of soil C0 2 is not necessarily the same as 
that of atmospheric C0 2 , although it has been assumed 
so in previous studies using 14 C to date soil carbonate. 
However, the present-day analytical precision at best is 
±3%o. and for most accelerator labs the analytical preci- 
sion is ±8%o. Therefore, when F is small, for example F 
= 0.1. the o l4 C value of soil C0 2 can be considered the 
same as that of atmospheric C0 2 . 

2) The 6 ,4 C values of soil CO : are not constant with depth. 

3) The 5 ,4 C values of soil C0 2 depend on the <5 14 C values 
of soil organic matter, soil respiration, and the relative 
contribution of C0 2 derived from organic matter decom- 
position to total C0 2 production ( F) . 

4) When the fraction of C0 2 derived from organic matter 
decomposition (F) * 0, the 6 14 C values of soil C0 2 vary 
continuously from the atmospheric value at the soil-at- 
mosphere interface to more negative values at depth. In- 
creasing the value of F, the 6 14 C values of soil C0 2 at 
any depth become more negative. 

5 ) When F = 0, the 6 ,4 C values of soil C0 2 vary continuously 
from the atmospheric value at the soil-atmosphere inter- 
face to more positive values at depth. This increase in 
5 ,4 C values is a result of diffusion effects on the different 
isotopic species (Cerling et al., 1991 ). If the t4 C data 
are corrected for isotope fractionation using l3 C data ( i.e., 
A ,4 C in Stuiver and Polach, 1977), this diffusion 
component will not be seen. 


Fig. I. 5 ,4 C value of soil C0 2 vs. depth for a model soil. It shows 
that the & l4 C values of soil CO: vary with depth from the ( pre-atomic 
bomb) atmosphenc value at the soil-atmosphere interface to more 
negative or more positive values at depth depending on the value 
of F. 


Figure 2 shows how the 5 14 C values of soil C0 2 vary with 
soil respiration rates. It is evident that at a given value of F, 
the higher the soil respiration, the more the <5 ,4 C values of 
soil C0 2 deviate from that of atmospheric C0 2 . 
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FlG . 2. Variations of the 6 !4 C values of soil C0 2 with soil respiration 
rates at a constant F value. The model shows that the higher the soil 
respiration rate, the more the 6 ,4 C values of soil C0 2 deviate from 
that of the atmospheric CO>. 


When all soil C0 2 is derived from root respiration and/ 
or decomposition of short-lived organic matter which has 
the same ,4 C content as the atmospheric C0 2 , the <5‘ 4 C values 
of soil CO: are relatively enriched (up to 8.5 %o at 15°C) 
compared to <5 14 C values of atmospheric C0 2 ( Fig. 3 ) due to 
diffusion effects. The lower the respiration rate or the higher 
the value of the diffusion coefficient, the less enriched the 
5 I4 C values of soil C0 2 are relative to the <5 l4 C values of 
atmospheric C0 2 . When soil respiration = 0, the <5 l4 C values 
Of soil CO: are the same as the 6 ,4 C values of atmospheric 
CO : . Again, considering the present-day analytical precision 
of ±8%o for I4 C analysis, this diffusion effect on l4 C content 
of soil CO : would be beyond detection. 

These figures show soil l4 CO: relationships for pre-indus- 
trial atmospheric conditions. To use this model to evaluate 
present-day soil conditions, one must change the atmospheric 
boundary conditions for C0 2 , I2 C0 2 , and l4 C0 2 . In our test 
of the model against empirical data given below, various 14 C aUT , 
values are used, based on the reported data. 

OBSERVATIONS IN SOILS 

It is important to establish if the above l4 C0 2 model is 
valid in soils. Presently, there are few 14 C0 2 depth profiles 
available. In Fig. 4, the l4 C data (reported as A l4 C, percent 
modem carbon (pmc) = (6 ,4 C/1000 + 1)(1 - 2*(25 
+ 6 ,3 C)/1000)* 100) for soil C0 2 and C0 2 concentrations 
(data obtained in May) from site #6 in Haas et al. (1983) 
and Thorstenson et al. ( 1983 ) are plotted (different symbols 
represent different sampling dates) and compared to our 
model calculation (solid line). Since there were no l4 C om . 
data and no C0 2 production information reported in their 
study, we assumed that the l4 C0 2 production is an exponen- 

f9C 14 

tial function of depth as <f^ 4 (r) = <^ 4 (0)e _bz with — — 

dz 

-Oatr = L(L - 380 cm was used in the model calculation), 
where <i>s 4 (0) is the production of l4 C0 2 at surface and b is 
a constant ( we use b - 0.003 in the model calculation ). We 


also assume that production of C0 2 is an exponential function 

SC* 

of depth in their soil as $f(z) = $?(())£ z(chz with — - 

dz 

= 0 at z = L' (L' = 520 cm was used in the model calculation ), 
where (0) is the production of C0 2 at surface and chz is 
the attenuation depth (we use c/iz = 300 cm in the model 
calculation). The solution to the diffusion Eqn. 2 with the 
assumed production function and boundary conditions are: 


C 4 = 


<*>s 14 (Q) 

(D\<)b 2 


(1 -e- zb )~ 


*i 4 (0)(z) _ 


(D“)b 


e-U’ + C" (9) 


<I> 5 *(0)(cAz 2 ) 
C ’ D • ( 


**> 


<*>,*(0 )(chz)(z) 

e L 

D* 


+ c?, do) 




FlG 3. & lA C value of soil C0 2 vs. depth for soils where all C0 2 is 
derived from root respiration and/or decomposition of short-lived 
organic matter. It can be seen that the 5 l4 C values of soil C0 2 are 
relatively enriched compared to 5 ,4 C values of atmospheric CO : due 
to diffusion effects. Figure 3(a) shows that the lower the respiration 
rate, the less enriched the 6 l4 C values of soil C0 2 are relative to the 
a‘ 4 C values of atmospheric C0 2 . Figure 3 (b) shows the effect of 
varying diffusion coefficient on the 14 C content of soil air. 
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Distribution of l4 C 


in soils 


where Cl 4 and Cf are l4 C0 2 and C0 2 concentrations in the 
soil; Cl* and C* are corresponding concentrations in the 
atmosphere; and L and L* are the depth of the non-flux 
boundaries for l4 C0 2 and C0 2 , respectively (the non-flux 
boundary is a boundary where C0 2 or ,4 C0 2 concentration 
gradient equals to zero). Other assumed parameters used in 
our model calculation are listed in Table 2. In their study, 
Thorstenson et aJ. ( 1983) concluded that the measured 
,4 C0 2 profiles cannot be readily explained with a diffusion 
model due to some as yet unexplained mechanisms. However, 
Fig. 4 suggests that our curve fits their data reasonably well, 
except for one datum point at the 15.9 meter depth. This 
point represents a groundwater C0 2 sample and, therefore, 
may be contaminated by carbon from other sources. The 
goodness of the fit is adequate to strongly suggest that vertical 
diffusion is indeed the dominant mass-transport mechanism 
affecting ,4 C0 2 and C0 2 distribution in this soil. It should 
be noted that 6 I4 C value of atmospheric C0 2 (the upper 
boundary condition ) is greater than 0%o (the boundary' con- 
dition used in Figs. 1-3). This is because present atmospheric 



Site#6 

Forest soil 

Pasture soil 

Respiration rate (mmoleVm2/hr) 

2 

n 

10 

Temperature ( B C) 

15 

27 

27 

Pressure (atm) 

1 

i 

1 

Porosity 

0.5 

0.55 

Cl 

P” 


0.66 

— m 

Diffusion coefficient (cm*/ sec) 

0.04 

0.095 

009 

« l3 Com 

-26%* 

-26%* 

-26%„ 

|4 C of respired CO? (PMC) 

129 



of atmospheric CO 2 (PMC) 

130 

114 5 

1 14 5 

Co? production function 

exponential 

constant 

constant 

Forest soil: A 14 C 0 m (percent modem carbon )=1 1 5 >0.005 ‘depth 
Pasture sod: A 14 ^,^ (percent modem carbon )= 100+0 035 ‘depth 
Atmospheric CO?: 350ppmv. 8 13 C=-7 


l4 C content is elevated above natural levels by nuclear weap- 
ons testing and use. 

Figure 5 compares the 14 CO : and C0 2 data from two Oxisol 
profiles (a forest soil and a 17-year-old pasture soil which is 
a degraded forest soil) from Paragominas in Brazil (Nepstad 
et al., unpubl. data) and our model calculations. The C0 2 
samples were collected in May, 1992 by Trumbore. Davidson. 
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Fig . 4. Comparisons of ( a ) C0 2 and ( b ) l4 C0 2 data from site #6 
m Haas et aJ. ( 1 983 ) with model calculation. Different symbols rep- 
resent different sampling dates and the curve is our model calculation. 
The goodness of the fits suggests that diffusion is indeed the dominant 
mechanism affecting the C0 2 and l4 C0 2 distribution in the soil. The 
datum point at 15.9 meter depth represents a groundwater sample 
and. therefore, may be contaminated by carbon from other sources. 



Soil C02 (%) 

Fig 5. Comparison of (a) ,4 C0 2 and (b) C0 2 (%) data from two 
Uxisol profiles from Paragominas in Brazil with model calculations 
The solid dots represent a forest soil and the open circles a pasture 
soil (a degraded forest soil). For each soil, the l4 C0 2 and C0 2 data 
are fitted using the same set of parameters ( i.e., same respiration rate 
diffusion coefficient, porosity, temperature, and pressure). Diffusion 
appears to be the dominant mechanism controlling the l4 C0 2 and 
C0 2 distributions in these two soils. 
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and Nepstad. The measured CO, respiration rates for May 
are 25 mmoles/m : /h for the forest soil and 10 mmoles/nr/h 
for the pasture soil. For each soil, the l4 CO, and C0 2 data 
are fitted using the same set of parameters ( i.e.. respiration 
rate, diffusion coefficient, porosity, temperature, and pres- 
sure) ( Table 2 ). The measured CO : respiration rate for the 
forest soil seems too high because it requires an unreasonably 
high diffusion coefficient and/or porosity to produce the ob- 
ser\ed soil CO : and CO : profiles. The high respiration rate 
reflects the decomposition of abundant litter on the forest 
floor which could have contributed significant amounts of 
C0 2 to the total CO ; flux measured at the soil surface and 
higher production of CO, in the upper 2 meters of the soil 
In our model calculation we used a value of 13 mmoles/ nr/h 
for respiration rate for the forest soil which gives a reasonable 
fit to both the M CO, and CO : data. Again, vertical diffusion 
appears to be the dominant mass-transport mechanism af- 
fecting IJ CO, and CO, distribution in these soils. However 
unlike the data of Haas et al. f 1983 ). the ,4 C content of soil 
CO : is greater than present atmospheric values. Since soil 
CO., is produced by root respiration and decomposition of 
soil organic matter which is a heterogeneous mixture of com- 
pounds turning over at different rates, the observed l4 C pro- 
files here suggest that the fractions of soil oreanic matter whose 
turnover contribute the most to the soil CO : contains a con- 
siderable quantity of bomb-produced IJ C. 

IMPLICATIONS FOR ,4 C DATING 
OF SOIL CARBONATE 


content of soil C0 2 can be depleted or enriched relative to 
that of atmospheric CO, depending on various factors. Soil 
carbonate formed during the early development of a soil, 
when l4 C content of soil organic matter is about the same as 
that of the atmospheric CO,, could have l4 C ages younger 
than the true age of the pedogensis. On the other hand, car- 
bonate formed later on in a soil, could have l4 C ages older 
than the true age of the carbonate precipitation. Elsewhere 
we (Amundson et al.. 1993 ) explore in detail the effects of 
our diffusion /reaction model on carbonate l4 C ages. 

CONCLUSIONS 

The distribution of l4 CO : in soils can be described bv a 
diffusion-reaction model. The l4 C content of soil CO, varies 
with depth depending on many factors: the l4 C content of 
soil organic matter, the relative contribution of root respi- 
ration to total CO, production, soil respiration rate, atmo- 
sphenc CO: concentration and ,4 C content, soil properties, 
temperature, etc. Our model suggests that l4 C ages of pedo- 
gemc carbonate could be older or younger than the true age 
of pedogensis. While there are other factors that can also 
affect the C age of a carbonate sample (such as sample 
thickness, i.e.. the total time required to form the sample 
being measured), our basic understanding of the initial l4 C 
contents of pedogenic carbonates is the first step to critically 
evaluating their potential as indicators of landform age. 

Editorial handling: E. J. Reardon 


Dates from soil carbonates have been considered unreliable 
estimates of the age of pedogensis because of unknown initial 
C/ C ratios in the carbonate and the possibility of subse- 
quent contamination with environmental M C (Callen et 
al.. 1983). Comparison of carbonate ,4 C ages with l4 C ages 
of coexisting organic matter suggests that radiocarbon dates 
calciflated from pedogenic carbonate in and areas were about 
-"00 to 7000 radiocarbon years too old ( Williams and Po- 
LaCH. 1969 ). On the other hand, radiocarbon dates of pedo- 
genic carbonate from the sub-humid pan of southeastern 
Australia were much younger than either the known age of 
deposition in which the carbonate is segregated, or the likelv 
age of pedogenesis ( Bowler and Polach. 1971). These dis- 
crepancies have been attributed to an initial low l4 C content 
of soil carbonate due to the limestone dilution effect and/or 
secondary contamination by environmental l4 C The lime 
stone dilution effect (Bartlett. 1951: Broecker and 
Walton, 1959) states that soil carbonate derives half of its 
C from dead calcium carbonate and another half from at- 
mospheric CO:, suggesting that radiocarbon age of such car- 
bonate would be about one half-life of 14 C (about 5570 years) 
older than the true age. However, studies on soil C0 2 and 
s C of soil C0 2 and pedogenic carbonates indicate that the 
CO: and its isotopic species in a soil system are in isotopic 
equilibrium. This implies that C derived from dissolving de- 
tntal carbonate is ultimately lost through isotopic exchange 
with soil CO, and isotopic composition of pedogenic car- 
bonate is determined by the isotopic composition of soil CO, 
(Cerling. 1984. 1991; Cerling et al.. 1989: Quade et al 
1989: Cerling and Quade. 1992). We have shown that l4 C 
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Carbon- 14 dates from soil carbonates have 
been considered unreliable estimates of the age 
of pedogenesis because of unknown initial 14 C/ 
12 C ratios in the carbonate and the possibility 
of subsequent contamination with environ- 
mental l4 C. Comparison of carbonate ,4 C ages 
with l4 C ages of coexisting organic matter sug- 
gests that radiocarbon dates calculated from 
pedogenic carbonate in arid areas were ~ 500- 
7000 radiocarbon years too old. On the other 
hand, radiocarbon dates of pedogenic carbon- 
ate from the sub-humid part of southeastern 
Australia were much younger than either the 
known age of deposition in which the carbon- 
ate was segregated, or the likely age of pedoge- 
nesis. These discrepancies have been attrib- 
uted to an initial low l4 C content of soil 
carbonate due to the “limestone dilution” ef- 
fect and/or secondary contamination by envi- 
ronmental l4 C. The “limestone dilution” ef- 
fect states that soil carbonate derives half of its 
C from “dead” calcium carbonate and another 
half from atmospheric C0 2 , suggesting that ra- 
diocarbon age of such carbonate would be 
about one half-life of l4 C (~5570 yr) older 
than the true age. However, studies on soil C0 2 
and <5 13 C of soil C0 2 and pedogenic carbon- 
ates indicate that the C0 2 and its isotopic spe- 
cies in a soil system are in isotopic equilib- 
rium. This implies that the ,4 C content of 
pedogenic carbonate should be determined by 


l4 C content of soil C0 2 and inherited “dead 
carbon” should not affect the age of soil pedo- 
genic carbonate. 

We present here a diffusion-reaction model 
that analyzes the l4 C0 2 distribution in soils. 
Carbon dioxide is produced in soils by biolog- 
ical processes and is transported to the atmo- 
sphere by diffusion. The diffusion mechanism 
applies to its isotopic species 1 2 C0 2 , 1 3 C0 2 and 
l4 C0 2 as well. Different isotopic species of C0 2 
react and diffuse independently of each other 
according their own concentration gradient and 
their own sources and sinks. Since the concen- 
tration of l4 C0 2 in a soil profile is controlled 
by the production and decay of l4 C0 2 , and by 
diffusion through the soil to the atmosphere, 
the concentration of l4 C0 2 can be described by 
a diffusion-reaction equation: 


dC' 4 d 2 C' 4 


6t 


= D\ 


di r 


+ 0i 4 -A C 


14 


(1) 


where C\* represents the l4 C0 2 concentration 
in the soil air ( mol cm ~ 3 ); D ] 4 is the diffusion 
coefficient of l4 C0 2 in the soil (cm 2 s - 1 ); 0* 4 
is the production of 14 C0 2 in the soil (mol 
cm -3 s _l ); and A the decay constant of l4 C 
(3.84* 10 -12 s"'). To model ,4 C0 2) informa- 
tion is needed regarding the production of 
l4 C0 2 (0i 4 ) in a soil. With the simplified as- 
sumption that soil C0 2 is primarily produced 
(a) by root respiration with practically no dif- 
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